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Abstract 

We project the Wilson/Pol chinski renormalization group equation onto its 
uniform external field dependent effective free energy and connected Green's 
functions. The result is a hierarchy of equations which admits a choice of 
"natural" truncation and closure schemes for nonperturbative approximate 
solution. In this way approximation schemes can be generated which avoid 
power series expansions in either fields or momenta. When following one clo- 
sure scheme the lowest order equation is the mean field approximation, while 
another closure scheme gives the "local potential approximation." Extension 
of these closure schemes to higher orders leads to interesting new questions 
regarding truncation schemes and the convergence of nonperturbative approx- 
imations. One scheme, based on a novel "momentum cluster decomposition" 
of the connected Green's functions, seems to offer new possibilities for accurate 
nonperturbative successive appproximation. 
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1 Introduction 



In recent years there has been a resurgence of interest in applying Wilson's exact renor- 
malization group (ERG) equation [U to problems in quantum field theory. Starting with 
Polchinski [[|, who derived Wilson's ERG equation in a context familiar to field theorists, 
it was first used to provide improved, simplified proofs of perturbative renormalizability 
for many field theoretic systems M. More recently, others have used this equation, or 
its Legendre transformed version 0, to study, in various nonperturbative approxima- 
tions, the qualitative, and in some cases quantitative, behaviors of many different systems, 
including scalar || |7j and vector || |[ field theories, gauge theories [JlQ|[ , fermionic sys- 
tems [11 1, and many phenomenological issues such as top quark and higgs masses fT2fl , 



bound states pi, and various aspects of QCD [fH|. Recent reviews may be found in 



15, 16, 17 



The challenge of nonperturbative approximations has been to reduce the infinite di- 
mensional space of couplings generated by the ERG equation to some computationally 
manageable size. This is often initially done simply by truncating the space to a few 
operators (usually the relevant and marginal operators of the canonical theory), and use- 
ful qualitative and even quantitative results are obtained by this approach. Improved 
approximations can then be sought by enlarging the space of effective interactions in 
some systematic way with the hope that results will converge towards the exact result. 
However, without a small parameter, as in perturbation theory, to give some measure 
of control over the degree of approximation, the degree of "improvement" is dificult to 
evaluate. In addition, additional nonphysical solutions are often generated as the finite 
space of couplings is enlarged [11], [L8|, [L9| , making it often difficult to identify with any a 
priori justification the physical solution sought. 

Another approach, which from the start includes an infinite number of interactions, 



is the derivative expansion |20| , [21], |22| , essentially a functional power series expansion of 
the Wilson effective action in powers of momenta coupling the Fourier transformed field 
variables. In this way all powers of the field are included at each level of approximation. 
Its lowest order of approximation, known as the "local potential approximation" (LPA) 0, 



23], [2^[ , is quite successful in providing qualitative and even reasonably good quantitative 
information regarding phase diagrams and critical exponents. It appears to generate 
no spurious solutions and permit the entire space of candidate renormalized theories to 



be searched for solutions in a simple, systematic way |15|, [19], |25fl. Its second order 
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of approximation gives impressive accuracy when compared with other approaches (e- 
expansion, exact, and high temperature expansion) |20, 21], and has offered the hope of 
leading to a practical method of successive approximation. 

A third approach, which has not been much used outside of perturbation theory, is 
complementary to the derivative expansion in that it involves an expansion of the Wilson 
effective action in powers of the field operators |]J, the expansion coefficients being the 
(zero-field) connected Green's functions of the theory. Thus, all powers of momenta are 
included at each stage of approximation. This approach has provided an alternative route 
to the e-expansion for scalar field theories |26] as well as, in the Legendre transformed 
version, a promising initial nonperturbative exploration of bound state questions [13] . 

In all of the above approaches the question of convergence of the approximation method 
involved is dificult to address a priori [||, [21], and although the approximations often 
improve when taken to their next higher order (usually the only one practical to compute, 
so far), there are indications that this good fortune cannot continue indefinitely. In 
|B| Morris studies the convergence of an expansion of the LPA in powers of the field 
variable. He demonstrates that the exact solution of the LPA equation has singularities 
in the complex plane that prevent the convergence of the power series expansion. If the 
presence of such singularities is a general feature of the ERG, one can expect it to prevent 
convergence of any expansion that involves a power series in the field operators. This 
would eliminate the hopes of convergence of the first and third methods discussed above. 
The derivative expansion also appears to have its own problem with convergence. In |p7f 
Dunne provides an exact solution for the effective potential for QED 2+1 in a particular 
inhomogeneous external magnetic field. From this he derives an all orders derivative 
expansion and shows that it is asymptotic. We take this to suggest that power series 
expansions in either fields or momenta are unable to provide a convergent approximation 
scheme for the ERG equation. 

It is interesting to contrast this situation with that of lattice real-space RG approx- 



imation schemes | 2q| , where a reliable (apparently convergent) method of successive ap- 
proximation seems to exist. One orders the space of interactions in terms of the interac- 
tion's relative degree of locality, i.e., nearest neighbor, next nearest neighbor, etc., and 
truncates accordingly. This locality property does not seem to be readily applicable to 
momentum-space formulations; it certainly has no relation to the power series approxima- 
tions discussed above, as all powers of both the field and momenta are typical components 
of local lattice interactions. 
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In this paper, motivated by the above considerations, we attempt to construct a 
method of successive approximation for Wilson/Polchinski's ERG equation for a scalar 
field theory in D dimensions that is not based on power series expansions. In particular, 
each order of approximation involves all powers of the field and (beyond the lowest order) 
all powers of momenta. In the final analysis it becomes essentially an expansion in the 
number of momenta coupling the Fourier transformed field variables, though it is, strictly 
speaking, in the space of connected Green's functions rather than the space of Wilson 
effective interactions that this expansion takes place. It is accomplished by projecting 
the ERG equation onto its uniform external field dependent effective free energy and con- 
nected Green's functions, creating an infinite hierarchy of partial differential equations 
(PDEs) to be truncated and solved numerically. A set of "fluctuation relations" (derived 
below) allows a choice of truncation schemes to be investigated. One of these schemes, 
based on a novel "momentum cluster decomposition" of the connected Green's functions, 
appears to have many advantages over the standard approach. Although it still remains 
to be numerically investigated, it appears to be a likely candidate for a convergent method 
of successive approximation as well as offering new possibilities for extending the power 
of the renormalization group approach in many directions. 

The remainder of this paper proceeds as follows. In Section 2 we introduce our system 
and notation and give a summary of the essential steps of the derivation of the Polchinski 
ERG equation and its Legendre transformed version. In Section 3 we derive a rescaled 
version of the Polchinski equation that is particularly well suited for the analysis that 
follows. We also show that it exactly reproduces Wilson's rescaled ERG equation. In 
Section 4 we derive a formalism for expressing the connected Green's functions in terms of 
the solution of the rescaled ERG equation. We demonstrate that our rescaling procedure 
gives the correct scaling laws for connected Green's functions. We also derive the set 
of "fluctuation relations" which will prove important in our later sections. In Section 
5 we derive flow equations for the connected Green's functions and present an analysis 
that suggests that the "obvious" truncation strategy will not be a convergent method of 
successive approximation. In Section 6 we present a new solution strategy based upon 
the "fluctuation relations" presented in Section 4. We define new "iV-point momentum 
clusters" and derive "momentum cluster flow equations" based upon a "momentum cluster 
decomposition" of the original connected Green's (iV-point) functions. We compare these 
new equations with the equations of Section 5 and discuss the apparent advantages of this 
new formulation. In Section 7 we offer some concluding comments. 
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One final note: Although in this paper we treat the calculation of the Gibbs free 
energy and connected Green's functions via the Wilson/Polchinski ERG equation, our 
entire approach can be equally applied to the calculation of the effective potential and one- 
particle- irreducible (1PI) vertex functions via the Legendre transformed ERG equation 
0, [5|]. To reflect this more general utility we use the more general terminology "AT-point 
functions" in the title and section headings of this paper. 



2 Flow Equations 

Our system is a self-interacting scalar field, <f)(x), coupled to an external source, J(x), 
in D Euclidean dimensions. We obtain the connected Green's functions from the IR- 
cutoff, t/V-regulated generating functional W^°[J] evaluated via the Wilson/Polchinski 
renormalization group equation. The derivation of this equation and of its Legendre 
transformed counterpart is, by now, fairly standard || [R|. For completeness, and to 



establish our notation, we will simply summarize, following Ellwanger flT3|] , the essential 
steps. 

The bare theory is regularized in the UV via a cutoff A in the propagator. We further 
introduce a running infrared cutoff, A, and corresponding propagator 

Pttf) se (^(g 2 ))- 1 = K ^ 2 ) - K ^ 2 ) (1) 

qz _|_ m z 

with 

K k (q 2 ) -> 1 for q 2 < A 2 

K A {q 2 ) -> for q 2 > A 2 . (2) 

The generating functional for cutoff connected Green's functions can be represented 
as the functional integral 

exp(-W£[J)) = jv<j> ex P {-i(0X<V) - S A ?[0] + (J,0)}, (3) 
where (J, 0) = / J q 4>- q , L = 72nW $ d D q, and (fi q and J q are the Fourier transforms of « 



x 



and J{x) respectively. Here the kinetic term of the bare action is represented in terms 
of the inverse propag ator Ri°{q 2 ), with the remaining part of the bare action denoted by 

The result of the functional integration in (|3|) may be formally represented by 

e -<°M = p |(j,p a a °j)X p -^ 



with 

We define an effective interaction Sj nt [0,A] by 



e -s iat [<t>,A) =f J$ > e -s&m m ( 6 ) 

We obtain Polchinski's flow equation for Si nt [(j),A] by differentiating equation (|B|) with 
respect to A: 

dS int 1 fdPt(q 2 )f 5 2 S int SStoSSi 



(7) 



OA 2J q dA V5<p q 5<p^ q 5<p q 5^ 
Our initial condition is Si nt [0,A o ] = S^[(f)]. 

The generating functional for cutoff connected Green's functions, W^°[J], is related 
to S int [(f>, A] through 

Wfr[J\ = S int [PtJ,A] - i(J,P A Ao J). (8) 

It satisfies the flow equation 

dw£° 1 r dRj°{g 2 ) f 5 2 Wf 5W£° 6W£° \ 
dA ~~2J q dA \8J q 5J- q ~ 5J q 5J- q J ' ^ 

The cutoff effective action, r^°, is given by the Legendre transform of W^°[J], 

ri°M = <°[J] + (J^), (io) 

where ip q = 5W£°[J]/5J_ q . If one further splits off a bare kinetic part of r A °[<^], 

r^] = ^,<V) + ftM, (ii) 

one obtains a flow equation for T^ [<p\ , the generator of 1PI vertex functions, of the form 

8fJ»_l f^Mf^^l- (12) 



OA 2J q OA I ^ y <fy 9 <fy>_ 9 , 
with initial value 

f£M = s&M • (is) 

It is interesting to note that unlike the flow equations for Si nt and T^ , which admit 
well posed initial value problems in terms of the bare action 5' in °[0], the flow equation for 
does not. This is because the substitution = P^°J used in eqs. (£|) and @ gives 
VFa °[J] = at P^ Q ° = 0. For this reason we use eqs. (0) and (||) rather than eq. @ for 
the calculation of W^°[J]. 
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3 Rescaling 

We recall that Wilson's original renormalization group program consists of two steps. The 
first step, performed by the derivation given above, integrates out the high momentum 
components of the field, generating an equivalent effective action with reduced cutoff. The 
second step rescales the momenta back to the original cutoff scale and rescales the fields 
so that the resulting equations are cutoff independent and have a fixed point solution 
corresponding to the scale-invariant, massless renormalized theory. Bell and Wilson have 
shown p9j that for "linear" renormalization group equations, such as Wilson/Polchinski 
above, the field rescaling must be appropriately carried out in terms of the scaling di- 
mension of the field or the resulting equations will not flow to a fixed point. The scale 
dependence of the cutoff function P A °(q 2 ) must also be chosen correctly. Finally, we also 
transform to dimensionless variables, which, while not really essential, is customary at 
this stage. 

We parametrize our effective momentum scale by 



t = log 



Ao 



(14) 



A ' 



and our dimensionless rescaled variables are 




(15) 



and 




(16) 



where dj and are the scaling dimensions of J(x) and <p(x) respectively, 



dj 



^(D + 2-r,), 



(17) 



and 



~(D-2 + V ), 



(18) 



where 77 is the anomalous dimension. A good choice for P A °(q 2 ) is 




(19) 



where 



P{q 2 /A 2 ) -> 1 for 
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Finally, by denning 

we get 



P(q 2 /A 2 ) -> for q 2 > A 2 . (20) 
SUtfA^StoMMMt)], (21) 



at 



^ ^nt - / </y 



l(D + 2-r])+q'-V g , 



50_ O ' 5(j) a , 5(j)_ g , J 



where 



(23) 

and the prime on V'/ means that it ignores the momentum conservation delta functions 
in 5S! mt /5<f)' ql . 

In the next section we will show that our choice of rescaling and cutoff functions gives 
the correct scaling laws at the critical point. Here we note that if we choose P{q 2 ) = e~ 2q2 , 
change variables to a' q = e q2 <p q , and define 

H[a',t) = -SU*V).*] ( 24 ) 

we get, after an integration by parts and neglecting all primes, 
dU r/D _ \5H 

I /(l 77 I 2g 2 ) fer — I — — I — ] (25) 
Jq 2 \ 9 5a q 5a q 5cr- q 5a q 6o~- q J ' 

which is exactly Wilson's equation jjfl. Eq. ( p2|) is more useful for our purposes, however, 
and we will continue to work with it in what follows. 



4 Scaling Laws and Fluctuation Relations 
for Appoint Functions 

In this section we derive a number of relations for connected Green's functions based on 
their representation in terms of the t — > oo (A — > 0) solution to the ERG eq. (|22). The 
[/^-regularized connected Green's functions are determined from the generating func- 
tional 

W A °[J\ = lim{s int [P A A V, A] - ^(J,P A A( V)}, (26) 



following eq. <^). Let 

00 1 r r 

Static*] = Ez? /,••• /, ^(g;,...,g;,t)0; r --0;,5 D (g; + --- + g;), (27) 
n=o n - J 1l q n 

where the above is written in terms of our dimensionless (primed) variables. Rewriting 
this in terms of our original (unprimed) variables, we get 

00 1 r r 

n=0 qi 1 n 

xA ~n djA -n( dj - d ,) K .. K5 D {qi + . . . + (2g) 

If we now substitute <f) q = P^°(q 2 )J q , we will then be able to study the limit A — > of 
eq. (PB|). However, because Po°{q 2 )J q = A.Q~ 2 P(q 2 / Al)J q is finite we can substitute it 
directly in ( p6|) to get 

00 -I 

W^[J] = limA^-/ •■■/ S' n (qiA-\... )qn k-\t)k- nd J 

A ^ u n=0 Hr! Jqi Jq n 

x P^/A 2 ) • ■■P(q 2 JAl)J qi ■ ■ ■ J q J D ( gi + --- + q n ) 

1 r Al 2 P(q 2 /Al)J q J_ q . (29) 



2 

In the same spirit we can now take the limit A — > 00, where A appears explicitly, to get 
the infinite cutoff generating functional 



00 1 r r 

W °°[J} = lim limA D ]T-/ •••/ ^(giA- 1 ,...,^- 1 ,*) 

A -»oo A^O ^ n! J 9 i J9„ 

xA- nd 'J qi ---J q J D (q 1 + --- + q n ), (30) 



where the remaining limit Ao — > 00 must be accompanied by appropriate fine tuning of 
the relevant bare couplings to produce an acceptable renormalized theory [[IJ [TjJ . 
The infinite cutoff connected Green's functions in a uniform external field j are 

G.fe, . . + ••• + «)=- ijT^W) (3D 

where Fq(j) = —Gq(J) is the Gibbs free energy, the Legendre transform of the effective 
potential. Thus, 

G k (qi,...,q k ,j) 

CO J 

= - lim limA^^- — S' n ( qi A-\...,q k A- l ,Q,...,0,t) 

A -+oo A-t-0 ^ (n — Kj! 

x A-k^A-*']"-*, (32) 
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where the G k and S' n are symmetric functions of their momentum arguments. If we define 

S k 



G k (q h ...,q k ,x,t)5 (qi H h q k ) 



TrS!Jcj>',t) 



<t>'=(27r) D xS D ( qi ) 



giving 



then 



G k (q 1 ,...,q k ,x,t) = - 



J k {n~k)\ 



S' n (qi,...,qk,0,---,0,t)x 



n—k 



G k (qi, q k ,j) = lim lim A D ' kdj G k (q 1 A" 1 , q k A~\ 3 A~ dj , t). 

Aq^oo A^O 



(33) 



(34) 



(35) 



This equation allows us to derive fixed-point scaling laws and asymptotic behaviors as 
follows. 

Let S{* t be a fixed-point solution of eq. (0): dS{* t /dt = 0. The corresponding G k are 

5 k 



G* k (q 1 ,...,q k ,x)5 (qi + --- + q k ) 



O:, 



9fc 



(36) 



In evaluating eq. (|3q ) for such fixed-point solutions, the absence of t-dependence allows 
the limit A — ► oo to be trivially taken, yielding 



Gt(qi, • • • , q k ,j) = lim A D -^^( ?1 A- X , . . . , g.A" 1 , jA~ 



A^O 



(37) 



For this equation to have a nontrivial limit we must have G k oc A~ D+kdj as A — > 0. 
We will consider this limit for two cases. First, for the case where ^=1,...^ 7^ 0, we let 
Xi = qiA~ l and y = jA~ dj . Then the Z{ = Xiy~ l l dj = qij~ l ^ dj are independent of A, and 



• x, 



- d J+17 a -D+kdj ~ d J J ~ 



A - gi - ...q k 
So for the limit A — > to exist we must have 

G* k (xi, ...,x k ,y) — ►&(*!,..., 



X 



r k 

k 1 



(3* 



(39) 



a: — kx) 
J/— ►oo 



where ^ is determined by solving the ERG equations. Then 
G* k (qi, • • • , q k ,j) = Gk(zi, • • • , z k )qi 



— —d + — 

k ■■■Qk 



(40) 



This is the general (hyper) scaling law for G k on the critical isotherm when 51=1,.. 7^ 



SKJ. For fc = 2, writing G 2 (q,-q,j) as G 2 (q,j), we have G* 2 (q,j) = G 2 (q/j D+2 -i)q 



77-2 



As this expression generally is used to define the anomalous dimension 77, its derivation 
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here validates our choice of rescaling operations and cutoff functions. For the second case, 
where j is finite and <&=!,...,& = 0, we must have 

G* k {0,...,0,y) — > constant • y~ h+ % . (41) 

Then 

G* k (0,..., 0, j) = constant • f k+ % . (42) 
From the definition of G k we also have 

n zeroes 

Gh(0i> ■ ■ ■ > ?fc-n> 0, . . . , 0, x, t) = — Gfc_ ft (gi, . . . , gfc-n, x, *), (43) 
which immediately gives 

. . . ,gjfe_ n ,0, . ..,0, j) = — G k _ n (q x , . . .,q k - n ,j)- (44) 

These relations will be important in what follows. They are a generalization of the so- 
called "fluctuation-dissipation theorem" or "linear response theorem" of statistical me- 
chanics, which relates fluctuations in a thermodynamic average to a corresponding suscep- 
tibility (linear response) [pOj , |3~l|j . The case most often encountered in statistical mechanics 



is for k = n = 2, which gives the well-known relation between the 2-point cumulant and 
order parameter susceptibility, 

G 2 (0,j) = d<p /dj, (45) 

where ipo — dGo/dj. We will refer to relations fl43 ) and (fPj) as fluctuation relations in 
what follows. 

5 Flow Equations for Appoint Functions 



The flow equations for the G k can now be derived by applying the definition (^J) to the 
flow equation (P2"|). Using ( |3"3"D we can write 

^i + --- + ^ = - L, ** • (46) 



9t V^---^ <9t =i2 7r)Dx6°( qi ) 

Then, by substituting the r.h.s. of eq. ( |22"1) for dS' int /dt and carrying out the indicated 
operations, we get (suppressing q, x, and t arguments when not necessary for clarity) 



dGn ~ , dGn . . . / dGn \ ^ /" , , . ~ , , 
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dG 



k>0 



at 



dG k k 

(D - kdj)G k - dj ~ 9i ' ^qfik 



i=i 



xG k - j+1 (-q ij+1 q ik , q ij+1 , . . . , q ik ) 

+ A{q')G k+2 {q',-ci 1 q 1 ,...,q k ), 
Jo 1 



-Qi4,Qii, •••»?*,) 



(47) 



where J 3 stands for the set of indices {zi, . . . , ij}, Z k stands for the set of consecutive 
integers {1, ...,&}, S/,-cx fc stands for the sum over the (*) different sets, Ij, of indices 
chosen from the set Z k , and the set {ij+i, ... ,4} is the complementary set Z k — Ij. We 
note that in the equation for dG k / dt we have qi + • ■ • + q k = from the 5-function of 
Using this together with eq. (fl"3| ) we can rewrite eqs. (f47j ) as 



5Go 

dG k>0 
dt 



~ dG / c)C \ 2 r ~ 

DG -djx-^+A(0)[^.j +j q A(q)G 2 (q) 

~ / dG \ dG k ~ 

(D - kdj)G k + (2A(0)-^ -djxj-^-^qi- V qi G k 

k-l , 

+ E £ 1 ^(fti + • • • + % (-?<i 

3=1 IjCZk 1 

xGi-j + i(-9i j+1 Qi k ,Qij + i, ■ ■ ■ , Qi k ) 



+ / A(g')G? fc+2 (g', -9', ft, 9fc ) 

Jo' 



(4f 



where we have explicitly pulled out the first and last terms from the double sum. 
The first four equations are 
dG 



dt 
dG 2 
dt 



dt 



dG 4 



~ f)G / dG \ 2 r ~ 



dG 



\dG 2 



(D - 2dj)G 2 + [2A(0)-g£ -djx^-q- V g G 2 
+ 2A(q)G 2 2 + [ A(q')G 4 (q',-q',q,-q) 



dG 



\dG 3 



(D - 3dj)G 3 + (2A(0)^ -djxj-^-Y^Qi- V 9l G 



r 3 



+ 2 



J2A( qi )G 2 ( qi ) 



i=i 



G 3 (qu 92, 93) + / A(q')G 5 (q', -q', q u q 2 , q 3 ) 

Jqi 



dG 



(D - Adj)G 4 + 2A(0)— - -djx 



dx 
11 



n , "\ dG 4 x - n 



1 



dx 



i=i 



r 4 



+ 2 



i=i 



£4(91,92,93,94) 



+ 2 ^(fti +fe) G '3(-9 J i - 9i2>9ii,9 J2 )G , 3(-9i 
^(9') ( 5e(9', ~9',9i, • • • ,94)- 



9i4 , 9«3 > 9m y 



(49) 



We note that Gi(q,x,t) = Gi(0,x,t) = dG (x,t)/dx due to the S D (q) in the definition of 
Gi(q, x, t), so we don't need to write a separate equation for it. Also, in response to the 
5-function with G 2 , we write (^2(9, —Q,x,t) as (^2(9, i). 

A natural procedure for solving the above hierarchy of equations is to truncate at 
a certain level by setting all higher order connected Green's functions equal to zero. 
To lowest order we get an equation for Go that has been studied long ago by Green, 



Gunton, and coworkers [[32|, |33|, |34| . Green |32] derived this equation via a steepest 
descent approximation to Wilson's functional integral form of his ERG equation ]l] . He 
showed that, with appropriate rescaling to ensure analyticity, mean field exponents are 
obtained f32"f . It was later shown that the equations reproduce more of the mean field 
phenomenology, including metastable and unstable branches in the free energy [34] and a 
spinodal fixed point p3 |. The higher order equations thus reflect systematic corrections 
to mean field behavior when successively higher order correlations are taken into account. 

A perspective on the higher order equations may be gained by realizing that they have 
already been studied for the zero-field case (x = 0) in D = 4 — e dimensions by Green 



and Shukla [25]. Because the equations are first-order quasi-linear they can be solved 
exactly in e-expansion, reproducing the standard results for the critical exponents. The 
e-expansion provides a natural truncation scheme for the hierarchy of equations, as Cm is 
of order e n_1 for n > 2. As the e-expansion is known to be asymptotic, we do not expect 
similar truncation schemes without the benefit of e as a small parameter to yield any 
improvement in convergence. Nor, due to continuity of the equations in x, do we expect 
different convergence properties for finite values of x. For this reason we do not expect 
the obvious truncation scheme to yield a convergent method of successive approximation. 
This is not to say, however, that low orders of approximation could not be useful. 



6 Momentum Cluster Flow Equations 

It is possible to improve upon the above truncation scheme by noting that, according to 
the fluctuation relations (|43|), each Gk contains derivatives of Gk><k within it. Thus the 
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truncation scheme discussed above needlessly throws away vital information. A truncation 
scheme that retains this information from higher order correlations can be constructed as 
follows. 

Due to the momentum conserving 5-function in the definition of G k , G k is a function of 
only (k — 1) independent momenta, and we will write it henceforth as G k (q\, ■ ■ ■ , q k -i, x, t), 
with the understanding that the omitted q k takes the value —q% — ■ ■ ■ — q k _ 1 . Defining the 
projection operator V gi G k (qi, ■ ■ ■ , qk-i, x, t) = G k {qi, ■ ■ ■ , qi = 0, . . . , q k -i, x, t), we define, 
for k > 2, 

r fe_1 

G k > 2 (qi,...,q k -i,x,t)= Y[(l-V qz ] 

L i=i 

We also define Gq(x, t) = Go(x, t) and G±(x, t) = 0. Clearly, G k >2 = if any of its (k — 1) 
momenta are zero. We call G k the k-point momentum cluster because all of its momentum 
variables are present in each term of its power series expansion. By expanding the product 
and making use of the fluctuation relations ([43D we obtain 



G k (qi,...,qk-i,x,t). (50) 



G k > 2 (qi, ■ ■ .,q k -i,x,t) 



k-2 gn _ f) k ~ 

EM)" E ^^_„,(g n ,...,g lfc _ 1 _„,x,t) + (-l) fc - 1 ^G . (51) 



n=0 ik-i 



In a similar manner, using the identity 1 = V Qi + (1 — V qi ), we find the inverse relation 
G k >2(qi, ■ ■ ■ , q k -i, x, t) 

k—2 gn ^ Qk ^ 

= E E ^Gk-n(q il ,...,qi k _ 1 _ n ,x,t) + —:G . (52) 

n=0 I k -l-n UX UX 

We call eq. (^) the momentum cluster decomposition of G k . 



Using these relations together with the flow equations for the G k (48) we can write 
flow equations for the G k as follows. Again we suppress unnecessary arguments for the 
sake of clarity. Define F k [G , . . . , G k ] by writing the flow equations ( ^8| ) as 



dG k 
dt 

Then 



F k [G ,...,G k ]+ [ A(q')G k+2 (q',-q', qi ,...,q k ^). (53) 

J a' 



9G n[ , liT d 2 G 



dt 



F [Go] + I a -qJ- + J A(q)G 2 (q) 
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dG 



k>2 



fc-1 



a 



dt 



n=0 



d 2 G k 
dx 2 



+ / A(q')G k+2 (q', -q', qi, ■ ■ ■ , qu-i) 



(54) 



where I a = J q A(q) and Fi[G ,Gi\ = dF [Go]/dx. To complete the derivation of the 
flow equations for Gk we must reexpress the arguments of F k _ n in terms of G k >< k - n via 
eq. (|52"D, but already in the above form we can see our most significant results, which we 
will describe shortly. 

The first four momentum cluster flow equations are 



dG 
dt 

dG 2 
dt 



dG, 
dt 



dG 4 
~dt 



DG — djx 



dGc 
dx 



+ A{0) 



dG \ 2 T d 2 G 



dx 



dx 2 



A{q)G 2 {q) 



(D-2 dj )G 2+ (2A(0)^ ' ~^ dG 



+ 2A(q)[G 2 + 



d 2 G 
dx 2 



(D - 3dj)G 3 + \2A{0) 

d^Gold^Go _ 
dx 2 [ dx 3 

G 2 (qi) + 



+ h 

dG 
dx 



dj x 

d 2 G 2 
dx 2 

- dj x 



) - 1 ■ V,G 2 - 2A(0) 
+ / A(q')G 4 (q',-q\q) 



d 2 G 
dx 2 



\dG 3 2 ~ 



+ 2A(0) 



+ 2 E M*) 



eC 2 



j^dG 2 ( gi ) 

i=i 
d 2 G n 



dx 



dx 2 



GziqxAi) + 



dG 2 ( gj ) dG 2 ( gi ) d 3 G 



dx 



dx 



dx 3 



+ 2A{q x + q 2 ) 



G 2 (qi + q 2 ) + 



d 2 G 
dx 2 



~ 2 dG 2 {q t ) d 3 G 
Gz{qi,q 2 ) + 2^ — ~ "+ 



i=i 



dx 



dx 3 



+ 1 



d 2 G, 



- I A(q')G 5 (q',-q',q 1 ,q 2 ) 
dGo 



dx 2 

(D - Adj)G A + (2A{Q) 
-d 2 G n ^ d 



— djx 



2A(0) 



dx 

E ^zG 3 ( qj ,q k ) 



\dG 4 3 ~ 
d 3 G * dG 2 ( gi ) , d 4 G d 2 G c 



dx 2 fri, dx 



E 



eC 3 



dx 3 dx 



<9x 4 dx 2 



+ 2 E Mod 



cC 3 



G 2 ( gi ) + 



d 2 G 
dx 2 



G 



q2, qs) 



d d d d 2 G (q ) d A G 

+ ^G 3 (q 3 , q k ) - -G z {q u Qj ) - -G 8 (ft, q k ) + + 
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d 2 G 2 (qA , d 3 G 



dx 2 



dx 3 



G 3 (qi,qj) + G 3 (qi,q k ) + 



3G 2 



dx 



+ 



dG 2 (q k 



dx 



o , , J n ( \ ^ gg^to) <9G 2 (gfc) <9 3 G 



E C 3 



X 



X 



G 3 (qj + q-k, Qi) + 



dG 2 



^G 3 ( qj ,q k ) + 



d 2 G 



dx 

z(qj] 



dx 



G 2 (qj + qk) + 



dx 



dx 3 



d 2 G 



dx 2 



+ 2A{ qi + q 2 + q 3 ) 



dx 2 
G 2 ( Ql 



d 2 G 2 {g k ) d 4 G 



q2 + q 3 ) 



dx 2 dx 4 

a 2 <5o 



{i,j,k} UX j=l 
eC 3 



dx 2 



+ 



dx 2 
d 4 G 



^4(91,52,93) 



&r 4 



d 2 G 4 
dx 2 



+ / A{q')G & {q',-q', qi,q 2 ,q 3 ), 



(55) 



where stands for the cyclic permutations of {1, ... , k}. 

The above equations display many features which appear to give them distinct ad- 
vantages over the standard Appoint flow equations (|49|). By comparing eqs. fl5"5j ) with 
eqs. fl49|) we see that the number of terms in the equation for dG k /dt increases dramati- 



cally with k relative to its "unclustered" counterpart (P9|). This is also directly apparent 
in comparing eqs. (]54]) and eqs. ([53]). The new terms are those which are discarded in the 
standard truncation approach to the flow equations (0), and their presence in fl5"3|) is one 
measure of the improvement of the momentum cluster approach. Because more terms 
are "rescued" as k increases this may offer hope for convergence of the approximation 
procedure. 

A more direct measure of the amount of information retained by the two sets of 
equations after truncation is obtained by comparing their integral terms, which contain the 
contribution from higher order functions which is lost when the equations are truncated. 
Their relative importance can be estimated by noting that, according to the G* k go 
to zero as an inverse power of the qi=i k as the qt go to infinity. As they generally are 
nonzero for qi = 0, we can assume their weight to be largest in the region near the origin 
and hence make a major contribution to the integral term in eqs. fl53|). This is because 
the function A(q) is typically concentrated near the origin as well, generally chosen to 
behave as e~ aq for large q. In contrast, the functions G k have had, by definition, their 
zero point values subtracted away, leaving functions that are concentrated away from the 
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origin, tending to nonzero values at infinity. These functions will therefore make a much 
smaller contribution to the integrals of eqs. (|4j) and, hopefully, be much less missed upon 
truncation. 

Another significant difference between the two sets of equations is that the flow 
equations for Gk are parabolic in x rather than first order, due to the presence of the 
lAd 2 Gk/dx 2 term in the equation for dGk/dt. This causes a dramatic change in their 
solution structure relative to the flow equations for the Gk- This is already apparent in 
that the new lowest order of approximation is the LPA approximation. This represents 
a significant improvement over its mean field counterpart of eqs. (|49|), both in terms of 
a more accurate calculation of critical exponents and description of universality classes 
and in terms of (presumably) avoiding the thermodynamically unstable branches of the 
mean field free energy. If the higher order equations of ([55]) represent as much of an 



improvement over their eq. (f|9]) counterparts as the LPA approximation is over the mean 
field approximation, we can expect these equations to prove very useful. 

Another consequence of the parabolic nature of the flow equations (|5SJ) is that they 
cannot be solved directly in the x = limit, as could eqs. (^). Thus, the full x-dependence 
of the flow equations will play an essential role in their approximation behavior. 

Of course, while the above points suggest the possibility of much improved approxi- 
mations based on truncations of the momentum cluster flow equations, the convergence of 
the truncation scheme remains an open question. Yet, a hopeful indication that conver- 
gence may not be out of reach lies in the fact that the momentum cluster flow equations 
are essentially an expansion of the ERG equation in terms of the number of momenta 
which couple the Fourier transformed field variables in the effective action. In a lattice 
real-space formulation this corresponds to an expansion of the effective lattice action in 
terms of local interactions ordered according to the number of interacting sites (of arbi- 
trary separation). Thus, we begin to make contact, perhaps with improvements, with the 
apparently convergent approximation scheme discussed in our Introduction. 

7 Concluding Discussion 

We close with some observations about future directions. The prospect for solving the 
lower order equations numerically looks quite good, as algorithms and approaches for 



solving PDEs are readily available p5[ and have already begun to be applied to the 



solution of approximate RG equations || [L3|, . These approaches should put solution 
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of the Gi equation within reach, at least in so far as it contributes to the equation for G2, 
as for that purpose only G&(qi, —qi, (fe, x, t) is required. Thus, at this level of truncation, 
only four variables, qf, q%, (gi + <?2) 2 , and x, are needed to solve for the fixed-point solution 
and critical exponents. To determine the full momentum dependence of G4, or the higher 
order corrections needed for the empirical study of convergence, will probably require 



further approximation [|13[ . 

Questions of convergence aside, the momentum cluster approach provides a possible 
advantage over the usual formulation in that it works directly with constant-field con- 
nected Green's functions rather than the interactions contributing to the RG effective 
action. Since the field is constant, it has no momentum index, considerably simplifying 
the description of momentum dependence when the approach is extended to vector models 
and gauge theories. This suggests that equations analogous to the G2 equation may be a 
useful starting approximation for the investigation of the space of renormalizable gauge 



theories, just as the LPA approximation was for scalar and vector models || ^5 



The momentum cluster approach can also be applied to the Legendre transformed 
ERG equation (|l"2D |4], |5|] and possibly to sharp-cutoff versions of the equation || [J7[ |38|] , 



as well. In this paper we limited ourselves to the Wilson/Polchinski version because it 
had the simplest form of nonlinear terms. More generally, one might consider whether 
the momentum cluster approach might have other applications, as many areas of physics 
make use of hierarchies of equations that begin with the mean field approximation and 
involve corrections from higher order correlations. If these equations can be represented 
in a generating functional format, the same fluctuation relations and momentum cluster 
decomposition might well apply, giving access to a new hierarchy of equations that no 
longer begins at a mean field level of approximation. 
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